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Israel-Stewart's causal theory of dissipative hydrodynamics, with the ADS/CFT lower limit of 
shear viscosity to entropy ratio (?7/s=0.08), give consistent description of a number of experimental 
observables in Au+Au collisions at RHIC (cm. energy ^/s=200 GeV) Assuming that in Pb+Pb 
collisions at LHC (cm. energy y / s=5.5 TeV), except for the initial temperature, other parameters of 
the fluid remain unchanged, we have predicted for the centrality dependence of multiplicity, mean pr, 
pr-spectra, elliptic flow. The central temperature of the fluid is adjusted to 7i=421 MeV such that 
in a Pb+Pb collision, with participant number N par t=350, average charge particle multiplicity is ~ 
900 and is consistent with the experimental trend observed at lower energies. Compare to Au+Au 
collisions at RHIC, in Pb+Pb collisions at LHC, on the average, particle multiplicity increases by 
a factor of ~1.6 , the mean pt is increased by ~10% only. The elliptic flow on the other hand 
decreases by ~15%. 

PACS numbers: 47.75+f, 25.75.-q, 25.75.Ld 
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I. INTRODUCTION 

Experiments in Au+Au collisions at RHIC 0, d, 0, 
5], produced convincing evidences that in non-central 
Au+Au collisions, a hot, dense, strongly interacting, col- 
lective QCD matter is created. Whether the matter can 
be characterized as the lattice QCD Q predicted Quark- 
Gluon-Plasma (QGP) or not, is still a question of debate. 
Relativistic hydrodynamics provides a convenient tool to 
analyse Au+Au collision data. It is assumed that in the 
collision a fireball is produced. Constituents of the fire- 
ball collide frequently to establish local thermal equilib- 
rium sufficiently fast and after a certain time Tj, hydrody- 
namics become applicable. If the macroscopic properties 
of the fluid e.g. energy density, pressure, velocity etc. are 
known at the equilibration time t,, the relativistic hydro- 
dynamic equations can be solved to give the space-time 
evolution of the fireball till a given freeze-out condition 
such that interactions between the constituents are too 
weak to continue the evolution. Using suitable algorithm 
(e.g. Cooper-Frye) information at the freeze-out can be 
converted into particle spectra and can be directly com- 
pared with experimental data. Thus, hydrodynamics, 
in an indirect way, can characterize the initial condition 
of the medium produced in heavy ion collisions. Hy- 
drodynamics equations are closed only with an equation 
of state and one can investigate the possibility of phase 
transition in the medium. 

A host of experimental data produced in Au+Au col- 
lisions at RHIC, at cm. energy -^1=200 GeV, have 
been successfully analysed using ideal hydrodynamics Q • 
Multiplicity, mean px, pr-spectra, elliptic flow etc. of 
identified particles, are well explained in the ideal hy- 
drodynamic model with QGP as the initial state. Ideal 
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hydrodynamics analysis of the RHIC data indicate that 
in central Au+Au collisions, at the equilibration time 
Tj w 0.6 fm, central energy density of the QGP fluid is 
It may be mentioned that ideal 
hydrodynamics description of data are not unblemished. 
Pt spectra or the elliptic flow are explained only up to 
transverse momenta ~ 1.5GeV. At higher px descrip- 
tion deteriorates. Also ideal hydrodynamic description to 
data gets poorer in peripheral collisions. 

Hydrodynamics implicitly assume local thermal equi- 
libration. In the rest frame of the fluid, particle mo- 
menta are isotropic. But at the early stage of the evo- 
lution, assumption of isotropic momentum distribution 
cannot be very accurate. System could only be in partial 
equilibration and for a creditable analysis of the exper- 
imental data, dissipative effects must be accounted for. 
Shear viscosity is the most important dissipative effect 
in heavy ion collisions. Shear viscosity of QGP matter 
is quite uncertain. In recent years, string theory moti- 
vated calculations § indicated that shear viscosity over 
entropy ratio of any matter is bounded from the lower 
side, rj/s > I/Att. It is then expected that at the mini- 
mum, shear viscosity over entropy of QGP matter should 
be 1/47T. 

In recent years there has been significant progress in 
numerical implementation of viscous dynamics [lj, 10, 
EH El, GS EI EE EE El El, El ■ At the Cyclotron 
Centre, Kolkata, we have developed a code " AZHYDRO- 
KOLKATA", to solve Israel-Stewart's 2nd order theory 
for dissipative hydrodynamics in 2+1 dimensions. In a 
recent publication 1], it was shown that minimally vis- 
cous hydrodynamics (rj/s— 0.08), with QGP as the ini- 
tial state consistently explain a large part of RHIC data, 
e.g. pT-spectra of identified particles, elliptic flow in 
minimum-bias/mid-central collisions etc. Indeed, the de- 
scription is better than that obtained in ideal dynamics, 
particularly at large pt- 

In not too distant future, Large Hadron Collider (LHC) 
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is expected to be operational. It is being planned to col- 
lide Lead beams at cm. energy ^=5.5 TeV, which is 
27 times larger than the RHIC energy collisions. One ex- 
pects that conclusive evidence for QGP formation can be 
obtained in Pb+Pb collisions at LHC. In the present pa- 
per, in the framework of minimally viscous hydrodynam- 
ics, we have given predictions for several experimental 
observables, which will be measured in Pb+Pb collisions 
at LHC energy. It will be shown that even if, from RHIC 
to LHC collisions, energy is increased by a factor of 27, 
the experimental observables e.g. multiplicity, pt spec- 
tra, mean px, elliptic flow etc. do show rather modest 
variation from RHIC energy collisions. The reason can 
be understood also. Multiplicity or entropy of the system 
increases logarithmically with energy. Thus, from RHIC 
to LHC, even though collision energy increases by a fac- 
tor 27, the entropy or multiplicity increase by a factor of 
~ 1.6 only. Experimental observables do not show very 
large variation from RHIC energy collisions. 
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centrality 


< Npart > 


< b > 


< Np ar t > 


<6> 


0-10 


346.8 


3.51 


318.9 


3.30 


10-20 


245.1 


6.13 


223.1 


5.79 


20-30 


169.9 


7.93 


154.1 


7.49 


30-40 


113.4 


9.38 


102.8 


8.87 


40-50 


71.6 


10.64 


65.0 


10.06 


50-60 


41.5 


11.76 


37.8 


11.11 


60-70 


21.5 


12.79 


19.8 


12.09 


70-80 


9.6 


13.76 


9.08 


13.0 


80-90 


3.1 


14.80 


2.97 


14.03 



TABLE I: Glauber model calculation for the average par- 
ticipant number < N par t > and average impact parame- 
ter < b > (in fm) for different ranges of centrality cuts in 
Pb+Pb/Au+Au collisions at LHC/RHIC energy. 



II. INITIAL CONDITIONS FOR PB+PB 
COLLISIONS AT LHC 

Details of our solution of Israel-Stewart's [2(| 2nd order 
theory of dissipative hydrodynamics can be found in [l[ . 
Briefly, assuming longitudinal boost-invariance, we have 
solved the Israel-Stewart's 2nd order theory for a baryon 
free fluid with dissipation due to shear viscosity only. In 
Israel- Stewart's theory, dissipative flows are treated as 
extended thermodynamic variables. For a baryon free 
fluid with only shear viscosity as the dissipative effect, 
energy-momentum conservation equation is required to 
be solved simultaneously with the relaxation equation for 
the shear stress tensor, 



= 0, (1) 
Dn^ = - — (n^ -2r 1 S7 < > 1 u u> ). (2) 

EqQ] is the conservation equation for the energy- 
momentum tensor, = (s + p)u^u v — pg^ v + , e, p 
and u being the energy density, pressure and fluid veloc- 
ity respectively. EqJ5] is the relaxation equation for the 
shear stress tensor tt^ . In Eq[2l D = u^d^ is the con- 
vective time derivative, V <fi u v> = i(Vu" + W) - 
^{d.u){g^ v — u^u v ) is a symmetric traceless tensor, rj 
is the shear viscosity and r w is the relaxation time. In 
[H Eqsffl and [2] are solved in (r = y/t 2 — z 2 , x, y, r] s — 
\ m £rf ) coordinates, assuming boost-invariance. 

We note that presently there is disagreement about 
the form of the relaxation equation to be used in heavy 
ion collisions. In [l5|, [l8| an extra term R = [u'V* + 
u v ir l,x ]Du\ is included in the relaxation equation. Shear 
stress tensor is traceless (ir 1 ^ = 0) and transverse to 4- 
velocity (u^ir^ = 0). The term is needed to maintain 
the transversality and tracelessness condition. Israel- 
Stewart [2(| developed the theory on gradient expan- 
sion of entropy, gradients of equilibrium thermodynami- 
cal variables are assumed to be small. The term [u tJ, ir l/X + 
u^ir^^Dux does not contribute to entropy and is missed 
in Israel-Stewart's theory. Moreover, in Israel-Stewart's 
theory, the term can be neglected (both and Du^ are 
small and their combination is neglected). We have also 
checked that for minimally viscous fluid, contribution of 
the term is negligible. For minimally viscous fluid, energy 
density evolution is hardly affected whether the term is 
present or not in the relaxation equation (see Fig. 8 and 
9 of Q). 

Details of the analysis of RHIC data in Au+Au col- 
lisions can be found in [l|. In brief, minimally viscous 
QGP fluid was initialised with a Glauber model initial 
condition, with 75% soft collisions and 25% hard colli- 
sions. In a b=0 collision, this corresponds to central en- 
tropy density £^=110 /m -3 , at the initial time Ti=0.6 
fm. The transverse fluid velocity at the initial time was 
assumed to be zero, v x = v y = 0. The shear stress ten- 
sor 7r Miy was assumed to attain boost-invariant value at 
the initial time T\. For the relaxation time tv, Boltz- 
mann approximation tv = 3?//2p is used. The freeze-out 
temperature was varied to fit elliptic flow in 16-23% cen- 
trality Au+Au collisions. It was seen that for Tp=130 
MeV, elliptic flow in 16-23% centrality Au+Au collisions 
as well as a host of other data are explained. For the 
equation of state we have used EOS-Q developed in 0], 
with bag model EOS for the QGP phase and hadronic 
resonance gas for the hadronic phase. EOS-Q has a 1st 
order phase transition at T c =164 MeV. 

We assume that in Pb+Pb collisions at LHC, except 
for the central entropy density, other parameters of the 
model remain unchanged. In Pb+Pb collisions also, the 
QGP fluid is thermalised at the same time as in Au+Au 
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FIG. 1: Filled circles are the PHENIX data for the charged 
particle multiplicity per participant -— ^ — -4p as a function 
of cm. energy for participant number N pa rt=350. The solid 
line is a fit to the PHENIX data by Eqj3] The unfilled circle is 
the extrapolated value of -g^J — at LHC energy y / s=5.5 
TeV, for participant number N par t = 350. 
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FIG. 2: The black, red and green solid and dashed lines are 
minimally viscous hydrodynamics predictions for the central- 
ity dependence of ir~ , K + and proton multiplicity in Pb+Pb 
collisions at LHC and in Au+Au collisions at RHIC. The 
yileds are normalised by a factor of N=1.4 to account for 
the neglect of resonance production. 



collisions at RHIC i.e. Ti=0.6 fm. The initial fluid veloc- 
ity is zero: v x = v y = 0, and the shear stress tensor has 
attained the boost-invariant value. And as in Au+Au 
collisions, in Pb+Pb collisions also, the hadronic fluid 
freezes-out at 7>=130 MeV. The central energy density 
or entropy density in Pb+Pb collisions at LHC energy 
(a/s=5.5 TeV) cannot be same as in Au+Au collisions 
at RHIC (^=200 GeV). One expects larger energy de- 
position in Pb+Pb collisions. To obtain the initial en- 
ergy/entropy density of the fireball in LHC energy colli- 
sions we proceed as follows: 

PHENIX collaboration [12] has tabulated the average 
charged particle multiplicity as a function of collision 
energy for a range of collision centrality. In Fig[TJ for 
participant number ./Vp Qrt =350, the average multiplicity 
_i — dN ah ■ snown as a function of collision energy . 
The multiplicity increases logarithmically with energy, 

^H=A + Bln^- S , (3) 
ai] 

with A — —0.33 and B = 0.75. We use the relation 
to extrapolate to LHC energy ^=5.5 TeV. The extrap- 
olated value of average charged particle multiplicity in 
LHC energy is ~ 927 + 70. We adjust the central entropy 
density to S ini —1%Q fm~ 3 such that a ^^=350 Pb+Pb 
collision produce ~ 900 charged particles. Entropy den- 
sity Si n i=180 fm~ 3 corresponds to central temperature 
T;=421 MeV. Compared to Au+Au collisions at RHIC 
(central temperature Ti=357 MeV), in Pb+Pb collisions 
at LHC, central temperature is ~ 20% higher. 

With the initial condition as described above we have 
solved the hydrodynamic equations and calculate invari- 
ant particle yield from the freeze-out surface at TV =130 
MeV. It may be noted that we are assuming boost- 
invariance. Consequently, our predictions are valid only 
in the mid-rapidity range. In the following, we will 
show our predictions as a function of collision central- 
ity or rather as a function of number of participants. 



Glauber model calculations for average participant num- 
ber < Np ar t > and average impact parameter < b > for 
different centrality cuts are given in table |H For Pb+Pb 
collisions at LHC energy, we have used <7« le i = 70 mb. For 
comparison, m table[H we have also shown the same re- 
sults for Au+Au collisions at RHIC, when ai ne i—4A mb. 



III. MULTIPLICITY AND MEAN p T IN PB+PB 
COLLISIONS AT LHC 

Particle multiplicity is one of the important observ- 
ables in heavy ion collisions. It is a measure of the en- 
tropy of the system. In Figf2l the black, red and green 
lines are minimally viscous hydrodynamic model predic- 
tions for the centrality dependence of 7r~, K + and proton 
multiplicity in Pb+Pb collisions at LHC. For comparison, 
minimally viscous hydrodynamic predictions for 7r~, K + 
and proton multiplicity in Au+Au collisions at RHIC are 
also shown (the dashed lines) in FigO We have neglected 
resonance production. To account for the neglect of res- 
onance production, yields are normalised by a factor of 
N = 1.4. From the predictions, it appears that compared 
to Au+Au collisions, in central/mid-central Pb+Pb colli- 
sions, particle yields are enhanced by a factor of ~ 1.6-1.8 
. It is expected. Multiplicity increases logarithmically 
with energy. 

Centrality dependence of mean px is another impor- 
tant observable. In Au+Au collisions at RHIC, in min- 
imally viscous hydrodyanmics, centrality dependence of 

< pt > of identified particles are reasonably well expa- 
lined. In Fig [3] minimally hydrodynamics predictions for 
mean px for 7r _ , K + and protons, in Pb+Pb colliisons 
at LHC are shown. In collisions beyond N part =50, mean 

< pt > is approximately constant; < pt >~ 0.5, 0.75 
and 1 for 7r~ , K + and protons respectively. For com- 
parison, predictions for mean px in Au+Au collisions at 
RHIC are shown in Fig[3]as the dashed lines. From RHIC 
to LHC, even though collision energy is increased by a 
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FIG. 3: The solid lines are predicted mean pt for ir~ , K + and 
proton in Pb+Pb collisions at LHC energy. For comparison, 
minimally viscous hydrodynamics predictions for mean pr in 
Au+Au collisions are also shown by the dashed lines. 



factor of 27, the mean px is increased marginally. For ex- 
ample, in a central collision with N par t=350, for all the 
species, mean pr is increased by ~ 10% from RHIC to 
LHC energy. 



IV. pt-SPECTRA IN PB+PB COLLISIONS AT 
LHC 

Minimally viscous hydrodynamic model predictions for 
7T- p T -spectra, in 0-10%, 10-20%, 20-30%, 30-40%, 40- 
50% and 50-60% centrality Pb+Pb collisions at LHC are 
shown in FigfJ] We have neglected resonance produc- 
tion. Resonance decay contribute to particle yield, more 
at low px than at large px- For example, at freeze-out 
temperature Tp=150 MeV, nearly ~ 50% of total pions 
are from resonance decay at Mt=0.5 GeV, contribution 
of decay pions decreases to ~20% at higher My=2 GeV 
[23l | . Minimally viscous hydrodynamics predictions, nor- 
malised by an factor of N=1.4, well explained the ir~ px 
spectra in Au+Au collisions. However, we must men- 
tion that overall normalisation do not account correctly 
for the resonance contribution to particle PT-spectra. pr 
spectra will be uncertain by 10-15%. For comparison, in 
FigHl predicted spectra in Au+Au collisions at RHIC are 
shown as the dashed lines, px spectra are slightly flat- 
tened at LHC. It is consistent with the predicted small 
increase in mean p? (see Fig|3]) in LHC energy. 

In FigJ5]and[5J we have shown the predicted px spectra 
for K + and protons. We have not shown, but here again, 
compared to RHIC energy pt spectra is flattened. Be- 
fore we digress, we would like to note that even though we 
have shown predictions right up to pt=5 GeV, at large 
Pt > 3 GeV, there may be other sources (e.g. jets) for 
particle production. It is unlikely that viscous dynamics 
will predict correctly particle production at high pt > 
3GeV. Present predictions are expected to be in reason- 
able agreement with future experiments up to px ~ 3 
GeV. 



FIG. 4: The solid lines are minimally viscous hydrodynamic 
predictions for tt~ p T spectra in 0-10%, 10-20%, 20-30%, 30- 
40%, 40-50% and 50-60% centrality Pb+Pb collisions at LHC. 
The dashed lines are the same for Au+Au collisions at RHIC. 
To account for neglect of resonance production, the yield are 
normalised by a factor of 1.4. 
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FIG. 5: The solid lines are minimally viscous hydrodynamics 
predictions for the K + invariant yield in Pb+Pb collisions at 
LHC energy. 



V. ELLIPTIC FLOW IN PB+PB COLLISIONS 
AT LHC 

One of the important observations in Au+Au collisions 
at RHIC is the significant elliptic flow in non-central col- 
lisions. Qualitatively, elliptic flow is explained in a hy- 
drodynamic model, re-scattering of secondaries generates 
pressure and drives the subsequent collective motion. In 
non-central collisions, the reaction zone is asymmetric 
(almond shaped), pressure gradient is large in one direc- 
tion and small in the other. The asymmetric pressure 
gradients generate the elliptic flow. As the fluid evolve 
and expands, asymmetry in the reaction zone decreases 
and comes a stage when the reaction zone become sym- 
metric and system no longer generate elliptic flow. Ellip- 
tic flow is an early time phenomena and a sensitive probe 
to the early stage of the fluid. 

In Fig[7J solid line is the predicted elliptic flow in mini- 
mum bias Pb+Pb collisions at LHC. For comparison, we 
have shown the minimum bias elliptic flow in Au+Au col- 
lisions at RHIC (the dahsed line). It is not shown here, 
but experimental minimum bias v 2 in Au+Au collisions 
is well reproduced in minimally viscous hydrodynamics. 
It is interesting to note that in Pb+Pb collisions, elliptic 
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FIG. 6: same as in Fig[5]but for protons. 
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FIG. 7: The solid line is the minimally viscous hydrodynam- 
ics prediction for the minimum bias elliptic flow in Pb+Pb 
collisions at LHC. For comparison, minimum bias elliptic flow 
in Au+Au collisions at RHIC is also shown (the dashed line). 



flow is reduced by ~ 15%. Reduction of elliptic flow in 
LHC has also been predicted by Krieg and Bleicher [25| . 
In a parton recombination model, they studied the energy 
dependence of elliptic flow in heavy ion collisions from 
AGS to LHC energy. It was observed that from RHIC 
to LHC energy elliptic flow decreases. Parton transport 
models also predict less elliptic flow in LHC energy than 
in Au+Au collisions at RHIC [11]. 

In Fig[SJ in four panels, we have shown the minimally 
viscous dynamics predictions for elliptic flow in 0-10%, 
10-20%, 20-30% and, 30-40% centrality Pb+Pb collisions 
at LHC. The black, red and green lines corresponds to 
elliptic flow for 7r~, K + and protons. Species dependence 
is similar to that at RHIC energy. At large pr, elliptic 
flow between different species are marginally different. 
Species dependence is seen only at low px, lighter the 
particle, more is the elliptic flow. We have not shown any 
comparison with vi at RHIC. But as with the minimum 



bias collisions, elliptic flow in different centrality ranges 
of collisions is reduced at LHC. 



VI. SUMMARY 

To summarise, in a minimally viscous (?7/s=0.08) hy- 
drodynamics , we have given predictions for several ex- 
perimental observables in Pb+Pb collisions at LHC en- 
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FIG. 8: (color online) Predicted elliptic flow in 0-10%, 10- 
20%, 20-30% and 30-40% Pb+Pb collisions at LHC are shown. 
The black, red and green lines are elliptic flow for tt~ , K + and 
protons respectively. 



ergy (y / s=5.5 TeV). Assuming that particle multiplic- 
ity increase logarithmically with cm. energy, we have 
extrapolated the lower energy data (tabulated by the 
PHENIX collaboration) to obtain an estimate of particle 
multiplicity at LHC energy collisions. Our estimate in- 
dicate that a Pb+Pb collision with participant number 
N par t=350, produces ~ 927 charged particles. The initial 
central entropy density of the fluid (Si n i— 180 fm~ 3 ) in 
Pb+Pb collisions was fixed to reproduce the extrapolated 
particle multiplicity. The initial time (r^) and freeze-out 
temperature (Tp) were kept fixed at the value, Ti=0.6 fm 
and Tp=130 MeV, as it was obtained in the analysis of 
Au+Au data in RHIC energy collisions. 

Our predictions indicate that compared to Au+Au col- 
lisions at RHIC, in Pb+Pb collisions at LHC, (i) particle 
multiplicity will increase by a factor of ~ 1.6, (ii) the 
mean p? will be enhanced by ~ 10%, (hi) pt spectra of 
identified particles will be (slightly) flattened and (iv) el- 
liptic flow will decrease by ~ 15%. We hope the results 
will be helpful in planning future experiments in LHC 
energy. 
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